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Recently, the CoGeNT collaboration presented a positive signal for an annual modulation in their 
data set. In light of the long standing annual modulation signal in DAMA/LIBRA, we analyze 

(T) the compatibility of both of these signal within the hypothesis of dark matter (DM) scattering on 

nuclei, taking into account existing experimental constraints. We consider the cases of elastic and 

i-C inelastic scattering with either spin-dependent or spin-independent coupling to nucleons. We allow 

for isospin violating interactions as well as for light mediators. We find that there is some tension 
between the size of the modulation signal and the time- integrated event excess in CoGeNT, making 
it difficult to explain both simultaneously. Moreover, within the wide range of DM interaction 
models considered, we do not find a simultaneous explanation of CoGeNT and DAMA/LIBRA 
compatible with constraints from other experiments. However, in certain cases part of the data 
can be made consistent. For example, the modulation signal from CoGeNT becomes consistent with 

^ the total rate and with limits from other DM searches at 90 %CL (but not with the DAMA/LIBRA 

^sO signal) if DM scattering is inelastic spin- independent with just the right couplings to protons and 

neutrons to reduce the scattering rate on xenon. Conversely the DAMA/LIBRA signal (but not 
CoGeNT) can be explained by spin-dependent inelastic DM scattering. 
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1. INTRODUCTION 



A number of experiments are searching with underground detectors for a signal from 
Dark Matter (DM) scattering on nuclei. A non-negligible interaction is for instance expected 
from DM that is a Weakly Interacting Massive Particle (WIMP). Two collaborations claim 
positive signals. Most recently, the CoGeNT experiment published results from 442 live 
days of data [1]. These confirm an event excess at low energies for the time integrated 
rate, already found in a previous release [2]. In addition, the data show a 2.8a indication 
in favor of an annual modulation of the event rate. In view of the long-standing evidence 
for annual modulation reported by the DAMA/LIBRA experiment [3] (DAMA for short) 
an intriguing question is, whether DM scattering can provide a consistent interpretation 
of the data from both experiments. The relevant WIMP mass region of m x ~ 10 GeV is 
also strongly constrained by other experiments, e.g., [4-6]. From a 4843 kg day exposure 
of the XENON100 detector, three candidate events have been observed, with a background 
expectation of 1.8 ± 0.6 events [6], showing no evidence for DM scattering. These data 
lead not only to the world strongest limit for conventional WIMP masses, constraining the 
DM-nucleon cross section to below 10~ 44 cm 2 for m x ~ 100 GeV, but have also important 
implications for the low mass region around 10 GeV. 

Motivated by these recent developments we give a critical discussion of the CoGeNT data, 
and reconsider possible DM interpretations of the global data, focusing on the hints coming 
from CoGeNT and DAMA. We entertain the possibility that one DM species is responsible 
for the two signals and vary the properties of the interactions with nuclei. We consider spin- 
dependent and spin-independent interactions for both elastic and inelastic scattering. We 
allow for general couplings to neutrons and protons, and we work in the contact interaction 
limit, but also analyze the effect of light mediators. Other analyses of the recent CoGeNT 
data have been performed in Refs. [7-10] and we compare our results to those studies in 
case of overlap. For related analyses prior to new CoGeNT data see, e.g., [11-26]. 

The outline of the paper is as follows. In sec. 2 we set the notation and define the DM 
interaction models used in the following. In sec. 3 we consider the simplest case, namely 
spin- independent elastic scattering and argue that there is tension between the modulation 
signal and the integrated rate in CoGeNT data, and we show that within this framework 
CoGeNT and DAMA are incompatible and disfavored by other constraints. In sees. 4 and 5 
we extend the analysis to inelastic scattering and generalized isospin dependence, for spin- 
independent and spin-dependent scattering, respectively. In sec. 6 we consider the nontrivial 
energy dependence introduced by light mediator particles. We conclude in sec. 7. 



2. NOTATION AND DEFINITIONS 



The differential cross section for non-relativistic DM scattering on a nucleon at rest, 

X N^ X N, is (N = n,p) 

do N 1 \M\ 2 
dq 2 647r m 2 N m x v 2 ' 

where v is the velocity of x an d the momentum transfer, q 2 , is related to the deposited 
energy Ed as q 2 = 2rriNEd. The nucleons are bound inside atomic nuclei with mass rriA, 
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in which case the momentum exchange is q 2 = 2rriAEd- In present-day DM experiments 
\q\ ~ (9(20 MeV) — 0(100 MeV). Assuming the same the corresponding \q\ for scattering 
on a single n or p would be \q\ ~ 0(0.2 MeV) - 0(1 MeV). 

In the paper we will consider also the case of light mediators active in the DM-nucleon 
interaction. The assumed masses will be in MeV range or even lighter. For concreteness let 
us consider the case of DM that is a fermion, interacting with visible matter by emitting a 
scalar of mass m^. The spin- averaged matrix element is given by 

2 2 

m x m N 



\M\ 2 = 1Q X 2 X X 2 N . * " (2) 

where X x ^n are the couplings of the scalar to %i N, 

Ant = X X <PXX + X n<PNN. (3) 

N=p,n 

The total cross section is obtained by integrating eq. (1) from q 2 = to the maximum 
value allowed by kinematics q 2 ^^ = with fi ab = m a m b / {m a + m b ). However, for a 

massless mediator the total cross section diverges. Therefore we introduce a minimal mo- 
mentum transfer qr^ in as an IR cut-off (corresponding to a threshold energy in an experiment) 
and define a reference cross section for DM-nucleon scattering as 

X ' ' 1 dq 2 n K + ^ im )K + 4/i> r 2 ef )' 



2 

min 



where we have neglected q^ in compared to q^ ax in the numerator. Motivated by typical 
kinematics in DM direct detection experiment, we choose q 2 ^ = (0.1 MeV) 2 and a reference 
velocity of v Te { = 10~ 3 • c. For /ijy ~ f^N (which holds for m x 3> m^r ~ 1 GeV) this 
corresponds to ~ (2 MeV) 2 . Note that in the limit of heavy mediators, m 2 3> 
the cross section becomes independent of g 2 ^ and t> re f, and is the total cross section for 
xN — >■ xiV scattering, as usual. 

The differential scattering rate in an experiment from x scattering off a nucleus (Z, A) in 
events/keV/kg/day is given by 

where p x is the local DM density and /(v) is the DM velocity distribution in the rest frame 
of the detector. The lower limit of the integration is set by the minimal velocity f mui that 
the incoming DM particle has to have in order to be able to deposit an energy Ed in the 
detector. For the case of inelastic xA —> x' 'A scattering it is given by 

^min = -=L== ( + o) , (6) 

where 5 = m x i — m x <C m x . The same equation also applies to elastic scattering, with 5 = 0. 
The above example of scalar interactions (3) leads to a spin independent (SI) scattering of 
DM on nuclei with the cross section given by 
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where q 2 = 2m aE^ and F p ^ n (q 2 ) are nuclear form factors. In general they can be different 
for neutrons and protons. In the following we adopt the approximation F n (q 2 ) = F p (q 2 ) = 
F(q 2 ), which corresponds to the assumption that neutrons and protons (and consequently 
charge and mass) have the same distribution within the nucleus. Experimentally, the charge 
distributions within nuclei are well known from elastic electron scattering. Recently, rea- 
sonably precise information on mass distributions for some nuclei were also obtained from 
coherent photoproduction of 7r° mesons. Within (9(5%) — 0(10%) experimental uncertainty 
the two distributions agree (see Table 3 of [27]), while there is some preference for smaller 
neutron root mean square (rms) radii. In contrast, mean field estimates predict larger neu- 
tron distributions with 0.05fm-0.2fm bigger rms for the heavy nuclei [28]. 

We use the Helm parameterization of the SI form factor, F(q 2 ) = 3e~ q s ' 2 [sin(/«r) — 
nr cos(nr)}/ (nr) 3 , with k = \q\ and s = 1 fm, r = \/R 2 — 5s 2 , R = 1.2A 1 / 3 fm [29]. We have 
checked that the Lewin-Smith parameterization of r, s [30] gives compatible results, with the 
difference for instance for XenonlOO bounds being less than 3% (see also the comparison in 
[31] of different form factor parameterizations and their effect on DM nucleon scattering). 

The coupling constants X p and A n are independent parameters of the model. If X p = X n 
the DM-nucleon interactions are isospin conserving. We allow for general isospin breaking 
interactions, A n ^ X p , and define 

tan# = ^, a = ^±^. (8) 
X p 2 y ; 

Further we have to take into account that for a given target Z several isotopes A* can be 
present. Neglecting the small effect of different isotopes on the kinematics and the form 
factors we define 

A 2 eS = J2 2r i [Zcos9 + (A i - Z)sin9] 2 , (9) 

iSisotopes 

where are the relative abundances of the isotopes. The SI cross section (7) can then be 
rewritten as 

'"^ A 2 eS G(q 2 )F 2 (q 2 ), (10) 



dE d 2/j, xp v 2 

where we used fi xp ~ \i xn . The function G(q 2 ) takes into account the effect of a light 



) = , ^ 2 • y iL ) 



mediator: 

[m^ + q 

In the limit of heavy mediators, ^> 100 MeV, this function goes to G(q 2 ) — > 1. If 
DM interactions are isospin conserving, X p = X n , the effective atomic weight is given by 
A 2 S — » A 2 = Yli r iA 2 , and one recovers the familiar expression for the SI cross section. 

The SD scattering is treated in an analogous way, allowing for isospin violating interac- 
tions. For scattering on a nucleus A with spin J we write 



dcrP m A 4ira SD 



dE d 2fx 2 xp v 2 3(2 J + 1) 



G{q 2 ) [a 2 S 00 (q 2 ) + ao^So^q 2 ) + a^iq 2 )} , (12) 



with 



tan 6=—, a = cos 9 + sin 9 , oi = cos 9 - sin 9 , a ST> = + a p B , (13) 



a p 
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a n ,a p are the couplings to the nucleons with af^ oc a^ p , and SV,- are nuclear form factors. 
Note the slightly different definition of a SD compared to the SI case in eq. (8), in order to 
obtain the conventional normalization in case of pure scattering off protons or neutrons. The 
SD form factors are computed according to ref. [32] for 133 Cs (abundant in the Csl crystals 
used by the KIMS experiment) and according to ref. [33] for all other nuclei. 

In the analysis we fix the astrophysics parameters and use everywhere a truncated 
Maxwellian velocity distribution in the galactic frame with escape velocity t> esc = 550 km/s, 
Vq = 230 km/s and the local DM density p x = 0.3GeV/cm 3 . The impact of astrophysical 
uncertainties was discussed for example in [15, 34, 35] for elastic spin-independent scattering 
and for the inelastic case in [36]. 

3. RECENT COGENT RESULTS AND ELASTIC SI SCATTERING 

We first discuss the recent CoGeNT results [1] in terms of the most common DM models 
- the ones leading to elastic scattering with spin-independent (SI) contact interactions, 
assuming f n = f p (this will be relaxed later). CoGeNT uses a germanium target with a very 
low threshold of 0.4 keVee (electron equivalent). To convert keVee into nuclear recoil energy 
measured in keV we use the quenching factor relation i^fkeVee] = 0.199(i?[keV]) 112 [37, 38]. 
Hence, the threshold of 0.4 keVee corresponds to a nuclear recoil energy of 1.9 keV, making 
CoGeNT especially sensitive to low mass WIMP scatterings. 

3.1. Modulation in CoGeNT? 

Let us evaluate the significance of the modulation signal reported in CoGeNT [1] , as well 
as whether a consistent interpretation in terms of DM together with the unmodulated rate 
can be obtained. 

Fig. 1 shows the CoGeNT data in the energy bins 0.5 — 0.9 keV, 0.9 — 3 keV, and 
3 — 4.5 keV. These data have been read off from fig. 4 of [1], where the 0.9 — 3 keV band has 
been obtained by subtracting the 0.5 — 0.9 keV data from the 0.5 — 3 keV data shown there. 
First we note that in all three cases the goodness-of-fit of no modulation is acceptable: we 
find x 2_va hies of 9.5, 16.8, 11.7 for 15 dof for the three bands (from low to high energies). 
However, if a cosine with a period of one year, an amplitude a, and a phase to is fitted, the 
X 2 for the two low energy bands improves significantly. We find 

0.5-0.9keV: A X l omod = 3.6, t = 117±29, a = 1.24 ±0.65, ( . 

0.9-3keV: &xl omo d = 9 -l > *o = 109 ± 18 , a = 0.56 ±0.18, ^ ' 

where to is given in day of the year and a in events/keV/day/kg. While for the lowest energy 
band the significance for the modulation is modest, at 1.9cr, for the 0.9 — 3 keV band we 
find slightly more than 3a (these values have been obtained by evaluating the A% 2 values 
for 1 dof). Note that the fitted to is not in perfect agreement with the expectation for a 
standard isotropic DM halo of to — 152 (June 2 nd ), especially for the 0.9 — 3 keV region, 
where it is about 2a off. 1 For the highest energy region 3 — 4.5 keV there is no significant 



A shift of the maximum of the DM scattering rate relative to June 2 nd can occur in more complicated 
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FIG. 1: Time binned CoGeNT data [1] with the mean values subtracted for three energy ranges. We 
show the results of fitting a cosine function with a period of one year and arbitrary phase (blue thin curve) , 
the best fit to modulation data only (green thick curve) with m x = 10.5 GeV, and the best fit to the 
unmodulated data (magenta dashed curve) with m x = 8 GeV. The corresponding x 2 values of the 0.5 — 0.9 
and 0.9 — 3 keV bands are also shown. 



preference for modulation. 

Now we proceed to a fit of these data in terms of elastic spin-independent DM scattering. 
We find a best fit point with m x = 10.5 GeV and o v = 9.5 x 10~ 41 cm 2 . The corresponding 
prediction for the modulation is shown in fig. 1 as a thick green curve. We observe the 
shift in the phase compared to the free cosine fit which leads to a worsening of the fit: the 
best fit DM predictions have only A% 2 = 2 (4.3) for 0.5 — 0.9 (0.9 — 3) keV compared to no 
modulation. The fit regions in the plane of DM mass and cross section are shown in fig. 2. 
The reduced Ax 2 values of the best fit point compared to no modulation imply that a closed 
region appears only at 90% CL, whereas at 99% CL data is compatible with no signal. 

In fig. 2 we show also the allowed region obtained from the unmodulated event rate. We 
use the data shown in the in-set of fig. 1 in [1], reproduced in the right panel. This is the 
unexplained event excess which remains after subtracting the L-shell peaks as well as a flat 
background. In the plot of [1] data points corresponding to two different peak-subtraction 
methods are shown (black and white data points). For our fit we take the average of the 
two points and conservatively use the lowest and highest edges of the error bars to account 
for this systematic uncertainty in the fit. For the region delimited by the magenta ellipses 
in fig. 2 we assume that the total event excess is explained by DM events. In this case the 
best fit point is obtained at m x = 8 GeV and a p = 7.4 x 10 -41 cm 2 , corresponding to the 



DM halos, see, e.g., [39]. 
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FIG. 2: Left: regions in the DM mass-cross section plane for elastic spin- independent scattering at 90% 
and 99% CL from CoGeNT data [1] using only the modulation data (green) or only the unmodulated rate 
(magenta). The corresponding best fit points are marked with triangles. For the fit to the unmodulated 
rate we assume that the total excess events (after subtraction of the L-shell peaks and a flat background) 
are explained by DM. The dashed curve corresponds to the 99% CL limit from CoGeNT unmodulated 
data requiring that not more events than observed are predicted. We also show the allowed regions from 
CoGeNT 2010 data [2]. Right: unmodulated CoGeNT event spectrum and two predictions from DM. The 
solid (dashed) curve refers to elastic (inelastic) scattering with parameters as given in the legend. 



solid curve in the right panel. The dashed curve in the left panel shows the limit which is 
obtained by just requiring that not more events are predicted from DM than observed. 

We see from fig. 2 that the best fit point from the modulation signal is excluded by 
the unmodulated rate. The allowed region from the unmodulated spectrum is still located 
within the 90% CL region of the modulated data, however, for those parameter values the 
modulation in the 0.9 — 3 keV region is significantly reduced. This is shown as dashed 
magenta curves in fig. 1. While the prediction in the lowest energy band is essentially 
identical to the one of the modulation best fit point, the modulation amplitude is reduced 
significantly for 0.9 — 3 keV due to the steeper shape of the energy spectrum at lower 
DM mass. The A% 2 compared to no modulation is now only 2.3, which means that the 
modulation is actually not explained in the region where the signal is strongest. 

A similar result is obtained for all parameter points to the left of the dashed curve in 
fig. 2. We conclude that in the simplest case of elastic spin-independent scattering within a 
standard halo there is tension between the modulation signal and the unmodulated spectrum 
in CoGeNT. 



3.2. CoGeNT vs DAMA vs CDMS vs XENON100 

Let us now discuss whether the CoGeNT results can be reconciled with other data in 
the framework of elastic spin-independent scattering. First we consider the DAMA/LIBRA 
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FIG. 3: Allowed regions (90% and 99% CL) and limits (90% CL) in the DM mass-cross section plane 
for elastic spin-independent scattering. We use data from CoGeNT [1] (modulated and unmodulated), 
DAM A/LIBRA [3] modulation, CDMS low-threshold Ge [4], CDMS Si [46], and XENON100 [6]. For 
DAM A/LIBRA we show the allowed regions for two illustrative assumptions on the quenching factor of 
sodium, for XENON100 we show the limit for two assumptions on the light-yield factor L e f[. 



data [3], which show annual modulation at significance of about 8.9<r. We fit the lowest 12 
bins in energy for the modulated event rate in DAMA, and use the unmodulated spectrum 
as an upper bound on the predicted rate, see [11, 34] for details. In the following we assume 
that ion- channeling does not occur, following the results of [40]. This leaves scattering on the 
Na atoms as possible explanation of the DAMA signal [41, 42]. Hence the quenching factor 
of sodium, g Na , is an important input into the analysis. In our fit we take into account the 
uncertainty of and marginalize the x 2 with respect to this parameters. In fig. 3 we show 
the DAMA allowed region for the two assumptions gNa = 0.3 ± 0.03 (the "default" value 
used in many analyses) and q^a = 0.5 ± 0.1. We have arbitrarily adopted such a high value, 
following the suggestion of [43] to reconcile CoGeNT and DAMA regions. We stress, however, 
that a recent measurement in [44] obtains a value of = 0.25±0.06 at E nr = 10 keV, which 
taken at face value strongly disfavors large quenching factors. In particular, the example of 
^Na = 0.5 ±0.1 used in fig. 3 is in clear conflict with the measurement from [44] and we adopt 
it only for illustrative purpose. Most of the older measurements compiled in Tab. 10 of [45] 
also seem to indicate a value gNa — 0.3. The DAMA allowed region is clearly separated from 
the region which can explain the CoGeNT unmodulated event excess, while there is some 
overlap of the region indicated by the CoGeNT modulation data. This holds even for the 
extreme assumption on the sodium quenching factor. 2 



2 Let us mention that our results differ here from the ones of Ref. [8] where a better agreement between 
CoGeNT and DAMA is obtained. The region from the CoGeNT unmodulated event rate in [8] is obtained 
at somewhat larger cross sections, whereas the DAMA region is found at cross sections of about a factor 



S 



Furthermore, fig. 3 shows that severe constraints from CDMS and XENON100 disfavor 
the region indicated by both, CoGeNT and DAMA. From CDMS we show in the figure 
limits from a dedicated low-threshold analysis of CDMS-II Ge data [4] using data down to a 
threshold of 2 keV, as well as data from a somewhat old data sample on Si [46] of 12 kg day 
exposure with a threshold at 7 keV. As discussed in [47], the energy scale determination is 
crucial for these limits. We have estimated that an ad-hoc shift in the energy scale of 40% 
may render the limit from [4] consistent with the CoGeNT region. Let us stress that we do 
not advocate such a shift, and the critique of [47] has been addressed in great detail in the 
appendix of [4] available in the arXiv version of that paper. 

XENON100 published results from a 4843 kg day exposure [6], obtaining three candidate 
events with a background expectation of 1.8 ± 0.6 events. We derive a conservative limit 
by using the maximum gap method from [48]. An important issue in the interpretation of 
data from xenon experiments is the scintillation light yield, conventionally parametrized by 
the function L eS (E nr ), see, e.g., [21]. An improved measurement of this function has been 
published recently [49]. For the XENON100 limit (thick curve) in fig. 3 we use for L e g (E nr ) 
the black solid curve from fig. 1 of [6], which has been obtained from a fit to various data, 
dominated by the measurement [49] at low energies. Such values of L e g are also supported by 
the analysis of Ref. [50] by modeling a relation between scintillation and ionization and using 
existing low energy data on ionization. In order to estimate the impact of the uncertainties 
related to the scintillation light yield we show also the limit (thin curve) which we obtain 
by using the lower edge of the light blue region in fig. 1 of [6] corresponding to the lower 
2a range for L e g. Nevertheless, we have also checked that the assumptions on L c r below 
E m m 3 keV (where no data are availble) have no impact on the XENON100 curve, and the 
limit remains unchanged even assuming L c g = below 3 keV. Hence, the result does not 
rely on any extrapolation into a region with no data. 

Fig. 3 shows that the XENON100 limit using the best-fit curve for L e fr excludes the 
CoGeNT (unmod) favored region, while pushing L e fj down by 2a leads to a region marginally 
consistent. In both cases large values for q-^ a have to be assumed to make DAMA consistent 
with the limits. Let us note that a dedicated analysis of XENON10 data [5] obtains the 
energy scale from the ionization signal and is therefore independent of the scintillation 
efficiency L e g. These results put additional severe constraints on the region of interest and 
exclude both the regions indicated by DAMA and CoGeNT (not shown). 

4. INELASTIC SCATTERING AND NONTRIVIAL ISOSPIN DEPENDENCE 

In the previous section we have seen that it is difficult to obtain a consistent interpretation 
of all data in terms of elastic spin- independent interactions. Here we investigate some 
possibilities beyond this simplest mechanism of DM interactions. 



two lower than ours. Our DAMA region agrees, e.g., with the one of Ref. [21]. 
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FIG. 4: Allowed regions and limits at 90% CL in the DM mass-cross section plane for inelastic spin- 
independent scattering with a DM mass spitting 6 — 25 keV. We use data from CoGeNT [1] (modulated 
and unmodulated), DAMA/LIBRA [3], CDMS low-threshold Ge [4], and XENON100 [6]. 



4.1. Inelastic scattering 

We start by considering inelastic scattering. This has been proposed in [51] in order to 
reconcile the DAMA modulation signal with the other constraints. For inelastic scatter- 
ing the modulation signal is enhanced compared to the unmodulated rate and in general 
heavy target nuclei are favored. While recent data from XENON100 [52] and CRESST-II 
(tungsten) [53] disfavor such an interpretation of DAMA, it might be interesting to use 
inelastic scattering to resolve the tension between the modulated and unmodulated rate in 
CoGeNT [7]. 

In fig. 4 we show the result for a DM mass splitting of 5 = 25 keV, which gives the best 
compatibility of modulated CoGeNT data with the unmodulated rate. As visible in the plot 
the best fit point for the modulation, m x = 12 GeV, a p = 2.2 x 10~ 39 cm 2 , is obtained just 
at the 90% CL limit of the unmodulated spectrum (used as an upper limit) as well as the 
CDMS Ge low-threshold bound. The predicted spectrum for the modulation as well as the 
X 2 is practically the same as for the elastic best fit point at m x = 10.5 GeV shown in fig. 1. 
Note, however, that inelastic scattering does not offer an explanation for the unmodulated 
event excess, as shown by the dashed curve in fig. 2 (right). The modified kinematics lead to 
a different shape of the spectrum which cannot explain the exponential event excess towards 
low energies. 

From fig. 4 we observe also that there is a clear mismatch between the DAMA and 
CoGeNT regions. Furthermore, the conflict with the XENON100 bounds persists, even 
in the case of using L e fj at the lower 2a limit. If the parameter space to the left of the 
XENON100 curve is considered, we find that again the modulation in the 0.9 — 3 keV 
CoGeNT data cannot be explained. We conclude that inelastic scattering does not provide 
a satisfactory description of all data. 
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FIG. 5: Left: A^ s /A 2 defined in eq. (9), as a function of the ratio of DM coupling to neutron and proton for 
various elements. Right: regions (90% and 99% CL) and limits (90% CL) in the DM mass-cross section plane 
for elastic spin-independent scattering with f n /f p = —0.7. The cross section on the vertical axis corresponds 
to a defined in eq. 8. We use data from CoGeNT [1] (unmodulated spectrum), DAM A/LIBRA [3], CDMS 
Ge low-threshold [4] and Si [46], XENON100 [6], and SIMPLE [56]. 



4.2. Generalized isospin dependence 

Motivated by the observation that Higgs-mediated SI scattering leads to roughly equal 
couplings of DM to neutron and proton (due to the dominance of the strange quark con- 
tribution), conventionally f n ~ f p is assumed, which leads to an A 2 factor for the SI cross 
section. This assumption needs not to be fulfilled in general. As pointed out recently in 
[24, 54] there might be negative interference between scattering off neutrons and protons, 
which would lead to a vanishing cross section for f n /f p = —Z/(A — Z) (for a realization of 
such isospin violating DM in a technicolor model see [55]). If a given element consists only 
of one isotope the cancellation can be complete, whereas if different isotopes are present 
only the contribution of one particular isotope can be cancelled exactly. 

Assuming natural isotopic abundances we show in the left panel of fig. 5 the suppression 
factor compared to the case f n = f p , where the effective atomic number in case of general 
couplings, Al s , is defined in eq. (9). We zoom into the region around the neutron/proton 
cancellation f n / f p ~ —1. One observes that due to the different neutron/proton ratios 
in the various elements, the cancellations occur at different values of f n j f p for different 
elements, which might allow better compatibility of some experiments. In particular, there 
is a minimum of the cross section for xenon at f n /f p ~ —0.7, which can be used to weaken 
the limit from XENON100 compared to CoGeNT (Ge) by more than one order of magnitude 
in the cross section [7, 24, 54]. 

Fig. 5 (right) shows the DAMA and CoGeNT allowed regions compared to limits from 
various other experiments for f n /f p = —0.7, chosen in order to weaken the XENON100 
constraint. Due to the enhancement of Na compared to Ge (see left panel) the DAMA region 
appears at lower cross sections and overlaps with the CoGeNT region at 99% CL, and indeed 
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FIG. 6: Allowed regions and limits at 90% CL in the DM mass-cross section plane for inelastic spin- 
independent scattering with f n /f p — —0.7 and DM mass splittings 6 = 15 keV (left) and <5 = 25 keV (right). 
The cross section on the vertical axis corresponds to a defined in eq. 8. We use data from CoGeNT [1] 
(modulation and limit from the unmodulated spectrum), DAMA/LIBRA [3], CDMS Ge low-threshold [4] 
and Si [46], XENON100 [6], and SIMPLE [56]. The best fit point of the CoGeNT modulation data is marked 
with a triangle. 

consistency with XENON100 is obtained, especially if the uncertainty on L e g is taken into 
account. Obviously, the CDMS low-threshold constraint from Ge cannot be circumvented 
due to the same material as in CoGeNT. Moreover, scattering on Si is enhanced compared 
to Ge (see left panel), which makes the conflict with the CDMS Si data more severe. In 
addition we show in this plot also a limit from the SIMPLE experiment [56]. This experiment 
uses a C 2 C1F 5 target and upper limits of 0.289 and 0.343 events/kg/day at 90% CL have 
been reported from 13.47 and 6.71 kg day exposures, respectively, with a threshold of 8 keV. 
The limits from SIMPLE are less stringent in the cases discussed so-far, but as visible from 
the left panel of fig. 5 the scattering on both CI and F is enhanced compared to Ge for 
f n /f p = -0.7, and in this case SIMPLE excludes the CoGeNT/DAMA region. 

We conclude that invoking a cancellation between neutrons and protons can lead to 
a consistent description of DAMA and CoGeNT (unmodulated excess), and weaken the 
XENON100 limit sufficiently However, the corresponding parameter region is excluded by 
CDMS Ge, Si and SIMPLE data. While we have shown results explicitly only for the choice 
fn/fp — —0.7 which minimizes the cross section for xenon, it is clear from fig. 5 (left) that 
also other choices for f n /f p cannot improve the global situation. 

4.3. Combining inelasticity and general isospin couplings 

It has been suggested in [7] to use inelastic scattering as well as f n / f P = —0.7 in order to 
reconcile the CoGeNT modulation with DAMA and other constraints. We reconsider this 
idea in fig. 6, showing regions and limits for DM mass splittings of 5 = 15 and 25 keV. As 
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FIG. 7: Allowed regions for DAMA/LIBRA [3] at 90%, 99% CL and limits from XENON100 [6], SIM- 
PLE [56], PICASSO [57], KIMS [58] at 90% CL in the DM mass-cross section plane for spin-dependent 
scattering off protons (tariff = in eq. 13). The left panel corresponds to elastic scattering, whereas the 
middle and the right panel correspond to inelastic scattering with DM mass splittings 6 = 22 kcV and 
S = 114 keV, respectively. 



discussed in sec. 4.1, for inelastic scattering the unmodulated spectrum in CoGeNT cannot 
be fit due to the specific spectral shape of the signal. Therefore we use the unmodulated 
rate only as an upper limit. For 5 = 15 keV the DAMA region includes the best fit point 
of the CoGeNT modulation, hence in this case the two modulation signals are in excellent 
agreement. However, the corresponding region is excluded by essentially all limits shown in 
the plot, including the one from CoGeNT (unmod). If the DM mass splitting is increased to 
25 keV, the CoGeNT modulation best fit point moves to smaller DM masses, becoming con- 
sistent with CDMS, SIMPLE, and CoGeNT (unmod). However, in this case the agreement 
between DAMA and CoGeNT (mod) is lost and the XENON100 limit becomes severe, since 
the heavy target is favored for inelastic scattering. Hence, a situation similar to inelastic 
scattering for f n = f p as shown in fig. 4 and discussed in sec. 4.1 is obtained. 



5. SPIN-DEPENDENT INTERACTIONS 



Let us briefly mention also the situation in the case of spin-dependent (SD) interactions. 
Here the relative importance of the various experiments is largely changed compared to the 
SI case, since the A 2 factor (or A 2 S in the case of general isospin dependence) is replaced 
by the corresponding nuclear structure function taking into account the spin-structure of 
the nuclei, see eq. (12). Depending on whether the spin of a nucleus is carried mainly by 
neutrons or protons there is a strong dependence on the isospin coupling, parametrized by 
tan# = a n /a p in eq. (13). 

The case of coupling only to the spin of the proton (9 = 0) is especially interesting 
for DAMA, since large part of the spin of Na and I are carried by protons. In this case, 
however, there are important bounds from experiments that use fluorine as a target: from 
COUPP [59], PICASSO [57], and most recently from SIMPLE [56]. As shown in the left 
panel of fig. 7 the fluorine bounds exclude the DAMA region. One possibility to make the 
DAMA explanation consistent with the bounds from the other experiments is to assume 
that the scattering is both spin dependent and inelastic, since then the rate for scattering 
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FIG. 8: Allowed regions (90% and 99% CL) and limits (90% CL) in the DM mass-cross section plane for 
elastic spin-dependent scattering with dominant coupling to neutrons (a n /a p sw 7). We use data from Co- 
GeNT [1] (unmodulated), DAMA/LIBRA [3] modulation, CDMS low-threshold Ge [4], and XENON100 [6]. 



on the light fluorine target is reduced. This possibility was discussed in terms of effective 
DM tensor interactions in Ref. [11] and in Ref. [19] by assuming that DM carries a large 
coupling to magnetic nuclear moments. As shown in the middle and right panels of fig. 7 
for DM mass splittings of 5 = 22 keV and 5 = 114 keV unconstrained allowed regions 
for DAMA appear. 3 These two solutions correspond to scattering off Na and I in DAMA, 
respectively. For large mass splittings the bound from the KIMS experiment (using a Csl 
target) become more relevant, but are still not able to exclude the DAMA region. Other 
constraints such as from Ge experiments or from the CRESST tungsten data disappear, 
since for those nuclei protons do not contribute to the spin. This, however, at the same 
time makes an explanation of CoGeNT data in this framework impossible. Constraints from 
neutrinos from DM annihilations in the sun have been considered for these DAMA solutions 
in [60]. 

For the germanium target in CoGeNT the spin is dominated by neutrons. However, in 
this case the constraints from DM scattering on xenon are particularly severe. By scanning 
the parameter 9 we have found best compatibility of CoGeNT and DAMA for 9 = 0.4457T, 
where a n /a p « 7, i.e., large coupling to neutrons. As visible in fig. 8, in this case the 99% CL 
regions of DAMA and CoGeNT (unmod) overlap. The corresponding region is excluded by 
CDMS low-threshold Ge data (which obviously follows the same pattern as CoGeNT) but 
also by the XENON100, due to the large contribution to the spin of the 129 Xe and 131 Xe 
isotopes from the neutron. 



3 Quantitative differences compared to the results of [11] appear due to slightly different values adopted for 
the halo parameters v and v csc . 
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FIG. 9: A\ 2 as a function of the mediator mass for various combinations of conflicting data sets within 
various DM interaction models. 

6. LIGHT MEDIATORS 

Up until now we have always assumed a contact interaction between DM and quarks, 
resulting from integrating out heavy mediator particle(s). "Heavy" means here masses large 
compared to the transferred momentum, i.e., m 2 ^> q 2 = 2m aE^ ~ (100 MeV) 2 . For masses 
smaller than this a nontrivial q 2 dependence will appear in the differential cross section. As 
a proxy we use the q 2 dependence that results from an exchange of light scalars, see eq. (7). 
In particular, for a "massless" mediator with m 2 ^ <C q 2 an additional factor 1/-EJ will appear, 
leading to an enhancement of the rate at low recoil energies. In this section we investigate 
whether light mediators could improve some of the fits discussed above due to the additional 
1/q 2 enhancement at low recoil energies. This could also have an additional benefit. Many of 
the considered models, such as inelastic scattering or neutron/proton cancellations, lead to 
a suppressed rate on the nucleus level, which has to be compensated by large DM-nucleon 
cross sections. One possible way to achieve such large cross sections is to consider light 
mediator particles. 

In fig. 9 we show the effect of decreasing the mediator mass for various combinations 
of conflicting data sets of the previous sections. For each combination we show the A% 2 
with respect to the x 2 minimum. Since we combine conflicting data sets (with the aim 
of improving the fit by reducing the mediator mass) the absolute value of the \ 2 is n °t of 
interest to us here. As is clear from the figure the fit does not improve for light mediators — on 
the contrary, in many cases the fit gets considerably worse. 

This is the case, for example, for the CoGeNT unmodulated excess versus DAM A (blue 
curve) and for DAMA excess versus XENON100 (red curve) assuming spin-independent DM 
scattering (cf. also fig. 3). Similarly, the conflict between the unmodulated CoGeNT rate 
versus the CDMS Si and SIMPLE constraints in the case of spin-independent scattering 
with f n / f p = —0.7 found in fig. 5 (black curve) becomes significantly worse for < 
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100 MeV. The same conclusion applies for the conflict between the CoGeNT rate versus the 
XENON100 bound in the case of spin-dependent interactions mainly on neutrons found in 
fig. 8 (orange curve). In all of these cases actually a lower bound on the mediator mass can 
be obtained, if the fit for large is accepted despite the tension between the data sets. The 
light mediator also does not remove the conflict between CoGeNT unmodulated spectrum 
versus CDMS low-threshold Ge data for SI scattering (cf. fig. 3), though in this case the 
X 2 remain basically flat compared to contact interaction. The same conclusion applies 
to the tension between CoGeNT modulation data and XENON100 bounds on inelastic 
scattering (cf. fig. 4), where again light mediators provide a fit of similar quality as the 
contact interaction. 

These results seem to be in disagreement with the ones of [10], where DM interacts 
via a massless mediator (in that case the photon) with the nuclei, and a consistent fit for 
CoGeNT, DAMA, and XENON100 is obtained. However, no direct comparison with our 
results is possible, since in Ref. [10] a multi-component DM model is considered with very 
peculiar halo properties due to strong self-interactions of the DM. 

7. CONCLUSIONS 

In conclusions, we have investigated the recent indication of annual modulation in the 
CoGeNT data, adopting the working hypothesis that this is a signal of dark matter scat- 
tering on nuclei. Assuming elastic spin- independent scattering, we find that there is tension 
between the modulation signal and the time integrated event excess in CoGeNT. For DM 
parameter values consistent with the unmodulated rate, the modulation amplitude in the 
energy range of 0.9-3 keV (where the significance of the modulation signal is highest) is 
too small to explain the data. Furthermore, we do not find any region of parameters where 
both CoGeNT and DAMA can be explained without being in conflict with bounds from 
XENON100 and CDMS. Motivated by these problems we considered a variety of more ex- 
otic DM interaction hypotheses. We tested inelastic spin-independent DM scattering, as 
well as elastic and inelastic spin-dependent DM scattering, allowing for isospin breaking in- 
teractions and light mediators. In all of the cases we find that it is not possible to reconcile 
both CoGeNT and DAMA at the same time, while also obeying existing bounds from other 
direct detection experiments. However, part of the data can be made consistent in some 
cases. For example, in a small parameter region the modulation signal in CoGeNT becomes 
consistent with the unmodulated rate (but cannot explain it, nor the DAMA signal), as well 
as with other bounds at 90% CL for inelastic spin-independent scattering with 5 ~ 25 keV 
and isospin breaking tuned to avoid xenon constraints [7] (cf. fig. 6). The DAMA signal 
(but not CoGeNT) can be due to DM if the scattering is both spin dependent and inelastic 
[11] (cf. fig. 7). 

As a general comment we mention that both the case of inelastic scattering as well as 
the case of cancellations between scatterings on neutrons and protons lead to a suppression 
of the event rate, which requires very large cross sections to account for possible signals in 
DAMA/CoGeNT, typically larger than 10~ 39 cm 2 . It remains to be shown that such large 
cross sections can be consistent with constraints from neutrinos [60, 61] and colliders [62-67]. 

Finally let us comment on preliminary results from the CRESST-II experiment, showing 
an unexplained event excess in their O-band data. They reported 32 observed events with an 
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expected background of 8.7 ± 1.4 [53]. A crude estimate indicates that such events could be 
consistent with elastic spin-independent scattering with similar parameter values as CoGeNT 
and DAMA [15]. For example, for m x = 10 GeV and a v = 10~ 40 cm 2 the predicted number 
of oxygen scatters for 400 kg day exposure with 100% efficiency is about 16 events, taking 
into account the reported individual thresholds of the 9 detectors. However, the indicated 
parameter region is again excluded by XENON100 and CDMS bounds [15]. For the model 
considered in fig. 6 we expect about 50 events for m x = 12 GeV, a p = 5 x 10 _37 cm 2 , and 
5 = 15 keV, whereas for 5 = 25 keV no scattering is possible on oxygen with CRESST-II 
thresholds due to kinematics. Similarly, also the spin-dependent inelastic scattering able to 
explain DAMA cannot accommodate the potential signal in CRESST, since 16 has no spin, 
and the inelasticity disfavors the scattering kinematically. 
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